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Abstract: In this study one resorts to the phenomenology of models endowed with a 
non-minimal coupling between matter and geometry, in order to develop a mechanism 
through which dynamics similar to that due to the presence of dark matter is generated. 
As a first attempt, one tries to account for the flattening of the galaxy rotation curves 
as an effect of the non-(covariant) conservation of the energy-momentum tensor of visible 
matter. Afterwards, one assumes instead that this non-minimal coupling modifies the 
scalar curvature in a way that can be interpreted as a dark matter component (albeit with 
negative pressure). It is concluded that it is possible to mimic known dark matter density 
profiles through an appropriate power-law coupling /2 = (i?/i?o)", with a negative index n 
— a fact that reflects the dominance of dark matter at large distances. The properties of 
the model are extensively discussed, and possible cosmological implications are addressed. 
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1. Introduction 



One of the main motivations of many modifications of gravity is the search for an alternative 
explanation for the problem of the flattening of the galaxy rotation curves [Q, ^, |3|, ^ . In 
the context of usual f{R) models 0], this has led to studies considering power-law 

curvature terms fi{R) oc in the action, instead of the linear curvature depicted in the 
General Relativity (GR) action. For n > 1, it is found that an addition to the Newtonian 
potential of the form A<^(r) = —GM{r/rc)^/2r, with /3 a function of n and being a 
parameter characteristic of each galaxy, may be used to match several rotation curves: the 
best fit indicates that f3 = 0.817, which arises from the choice n = 3.5 |8|. 

In another study, it is found that a metric describing a flat asymptotic orbital velocity v 
may be obtained by assuming that a logarithmic factor affects the curvature term, fi{R) = 
foR{l + f logi?); this may be approached by a power-law with an exponent close to unity 
(given the non-relativistic condition u <C 1), fi{R) ~ foR^^'^'^ Notice that this approach 
exhibits no galaxy-dependent parameters, and yields a universal asymptotic velocity v, in 
contrast with the Tully-Fisher and Faber-Jackson relations, which roughly-speaking posit 
a dependence of the visible mass M of a galaxy with a power v"^, for spiral (m = 4) and 
elliptical galaxies (m = 6), respectively. 

Usual f{R) models are mostly considered as low-energy phenomenological models that 
attempt to grasp some of the possible implications of a more complete theory: indeed, it is 
known that one-loop renormalization of GR requires the introduction of higher order terms 
in the curvature in the Einstein-Hilbert action functional; furthermore, other available 
invariants — such as contractions of the Ricci or of the Riemann tensor — may also arise 
when quantum corrections arising from string theory are considered (see Ref. [ p!o| ] for a 
thorough discussion). 

This said, it is interesting to also consider another fundamental aspect common to some 
extensions of GR: the existence of a non-minimal coupling between matter and geometry. 
Indeed, while GR presents us the minimal coupling ^—gCm (a tensorial density of weight 
1) in the action, more evolved forms may arise in the context of scalar-tensor theories (see 
discussion of Appendix B) if matter is represented by a scalar field, or due to the effect of 
one-loop vacuum-polarization in the formulation of Quantum Electrodynamics in a curved 
spacetime 

The structure of this study is as follows: in section |2|, the non-minimal gravitational 
coupling model [12| is introduced and its main features are discussed; session ^ presents 
a proposal to account for the flattening of the galaxy rotation curves, resorting to non- 
geodesical motion present in the model. In section ^ a more evolved scenario is studied, 
where the non-minimal coupling of geometry with a null dust matter distribution originates 
a dynamical effect that can be interpreted as dark matter. The dynamical features of this 
approach are then explored and, in section ^ applied to two well-known dark matter 
distributions: the Navarro-Frenk- White (NFW) cusped density profile and the isothermal 
sphere model. 
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2. The model 



Following the introductory discussion, one postulates the following action for the theory 



s 



h{R) + [l + \f2{R)]C„ 



-gd^x, 



(2.1) 



where fi{R) (with i = 1,2) are arbitrary functions of the scalar curvature R, Cm is the 
Lagrangian density of matter and g is the metric determinant. The contribution of the non- 
minimal coupling of /2 is gauged through the coupling constant A, which has dimensions 
[A] = [/2]~^- The standard Einstein-Hilbert action is recovered by taking /2 = and 
fi = 2k{R — 2A), where k = c^/IGttG and A is the cosmological constant (from now on, 
one works in a unit system where c = 1). 

Variation with respect to the metric g^i, yields the modified Einstein equations of 
motion, here arranged as 



(Fi + 2XF2Cm) R^u - ^flg^^u = {^^^^ - 9/.^n) [Fi + 2AF2£„,) + (1 + A/2) T^^, (2.2) 

where one defines O/iu = V^Vi, for convenience, as well as Fi{R) = f[{R), and omitted the 
argument. The matter energy-momentum tensor is, as usual, defined as 



2 5 (i/^£m) 

^ Sg 



111/ 



so that the trace of Eq. (|2.2D reads 



(2.3) 



(Fi + 2AF2£„) R-2h= - 3n (Fi + 2AF2£„) + (1 + A/2) T. (2.4) 
The Bianchi identities, V^Gf^u = imply the non-(covariant) conservation law 



AF2 



1 + A/2 

which is interpreted as due to an energy- momentum exchange between matter and geometry 
due to the non-trivial fi{R) and /2(-R) terms [ p^ . 

Since a full study of the joint effect of a non-trivial fi{R) and f2{R) is too involved, we 
focus our attention on the latter, thus setting fi{R) = 2kR (discarding the cosmological 



constant A); this reduces Eq. (2.2) to 



1 + -F2Cm) R^u - ^Rg^u = - (nfiu - g^un) {F2Cm) + 7^ (1 + A/2) T^u, 

Ki / A Ki AK 



and, from the trace, the equivalent of Eq. (p.4| 



(2.6) 



l--F2£„)i? = 3-n(F2£, 

K J K 



2k 



:i+A/2)T. 



(2.7) 



-3- 



3. Flattening of the galaxy rotation curves due to non-geodesical motion 

As discussed in section ||, the flattening of the rotation curves of galaxies has been studied 
in the context of a non-trivial curvature term fi{R) oc R"' These approaches share 

with several other similar studies the assumption that motion of test particles is geodesic. 
Clearly, although most extensions of GR "inherit" its metricity, which states that spacetime 
is endowed with a metric Qf^^, describing geodesies, motion is geodesical only if the energy- 



momentum tensor Tfiv is conserved — and, as Eq. (2.5) shows, this may not be always 
the case. Hence, one may first attempt to solve the problem of the flattening of the galaxy 
rotation curves through the effect of this non-geodesical nature of motion. 

One begins by assuming that matter is a perfect fluid, with an energy-momentum 
tensor described by 

T^u = {p + p)U^Uu+pgf,u, (3.1) 

where p is the energy density, p{p) is the pressure (given by a suitable equation of state) 
and is the four- velocity, satisfying the normalization condition U^JJ^ = —1; Eq. ( p.5| ) 
then translates into an additional force, 



P + P 



{Cra+p)VuR + yvP 



l + A/2 



h^'^, (3.2) 



with the projection operator hfj_y = g^y — U^Uy (so that h^^U^ = and the extra force is 
orthogonal to the four- velocity) . 

Classically, the Lagrangian density Cm of a perfect fluid has several equivalent, on-shell 
expressions: for instance, £mO = —p, ^mi = P and Cm2 = —na, where n is the particle 
number density, and a is the physical free energy, defined as a = p/n — Ts, with T the fluid 
temperature and s the entropy per particle. However, it was shown that, in the presence 
of a non-minimal gravitational coupling with matter, this degeneracy is lifted and one has 



to consider the original, bare Lagrangian density Cm = —p 114 1. 

Furthermore, one may consider the non-relativistic case, whereas matter is described 
by a pressureless dust distribution for which p <^ pc?; from a more rigorous viewpoint, the 
thermodynamical formulation of an action functional for a perfect fluid and subsequent 
relation p = n{dp/dn) — p implies that a null dust distribution is characterized by an 
equation of state of the form p (x n. The energy- momentum tensor simplifies to T^i^ = 
pUpUu, with trace T = —p. 

For simplicity, one takes a static, spherically or cylindrically symmetric case. This 
Ansatz and the above choices yield 

i + A/2 

where the prime denotes differentiation with respect to the radial coordinate r. 

Finally, one assumes that the velocity only has only a tangential component, Uq = v{r), 
with the time component Uq given by the normalization condition U^U^ = — 1 — which is 
valid for the spherical or cylindrical scenarios; this yields /i^'" = 5^^. 
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Before proceeding, one points out that the goal of this section suggests that the effect 
of the non- minimal gravitational coupling to matter is non-perturbative, Xf2{R) S> 1; this 
would yield a radically different set of equations of motion (Eq. (|2.6D ) than those of GR, and 
would prevent the classical identification R = —SttGT. However, this stems from a liberal 
use of what "perturbative" means. Indeed, one may reasonably aim for a non-perturbative 
modification of the geodesical motion — which depend on the factor XF2/{1 + A/2), while 
still maintaining a perturbative regime for the metric, that is, equations of motion which 
closely resemble those of GR. This is achieved if 



|A/2(i?)|«l , X[F2{R)p]" p, (3.4) 

having used T = Cm = —p- 

These two relations guarantee that Eq. (2/7) reads R ~ —SttGT = SirGp; assuming 
this, one may consider that the geodesical motion of a test particle is described by the usual 
Newtonian potential — given that the non-relativistic regime is valid for a typical galaxy, 
otv = GM/r^, M being the visible mass of the galaxy, assuming one is sufficiently far away 
from its center — while the non-geodesical force per unit mass/*" leads to the fiattening 
of the rotation curve; sufficiently far away from the center of the matter distribution, the 
latter should dominate, |/^| ^ a^- Assuming an asymptotic value v{r) — )■ Voo, one gets 

GM GM , , 

r >aiv^ — >^^r>^, 3.5 

which, using typical values M ~ 10-'^'^ Mq and v^o ~ 200 km/s, yields r > 1 kpc. 

In the Newtonian limit, valid for the region under consideration r > 1 kpc, one assumes 
a diagonal metric, and Eq. (|3.3| ) for an extra force that is purely radial. As previously 



discussed, one assumes that the field equations Eq. (2^) are only weakly perturbed, so 
that one may substitute R = SirGp; this, together with condition A/2 <^ 1 and g^"^ ~ 1 
yields 

r -XF^R'ir) = -XfMr)) = - — , (3.6) 

r 

after using the chain rule for differentiation; given that in the subsequent sections the 
possibility of non-horizontal fiattening (that is, non-vanishing slope for v'^{r)) will also be 
addressed, one considers instead a more general equation, 

\fMr)) = '^{-)\ (3.7) 

where r^, is an unknown normalization constant and j3 is the outer slope of the rotation 
curve. 

The solution then reads 



A/2(p(r)) = <^ 



vl^logU) ,/3 = 



(3.8) 
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where r* is used in the /3 = case as an integration constant ( i.e. its value is not shared 
between the two branches). Without loss of generality, one may set A = « 10~^^. 

In order to assess the dependence of /2 on the curvature, one simply needs the expres- 
sion for p{r) with the above equation; assuming that, sufficiently far away from the galactic 
core, the visible matter density profile takes the form p(r) ~ pQ{a/ry^ (with m ^ 0), one 
gets r = a a/po/p- Substituting into Eq. (3^), one has 



f2{p) 



.1 JL 

13 U. 



/3/m 



i-log(f) /3 = 



(3.9) 



after defining = po{a/r^)"^. Recalling that R ^ SirGp, one gets 



^ ( ,/3/0 



-i-log(^) ,/3 = 



f2{R) = { 

which exhibits the characteristic curvature = SvrGp*. This yields the condition 



(3.10) 



F2{R) 



1 ( R 



mR \R 



fi/m 



(3.11) 



for all /3. 

One must now show that the non-minimal coupling given by Eq. ( 3.1C| ) satisfies the 



perturbative conditions Eqs. (3.4): 



f (7) «1 ' /3^0 



log ( — 

r 



< 1 , /3 = 0, 



(3.12) 



and 



^^/3(l + /3) - - - «1, 



(3.13) 



with the additional lengthscale rg = 1/ \JSttGpq oc (where the latter is the Schwarzschild 
radius of the galaxy). For a given visible matter density profile p{r) characterized by values 
of ro and a, the former conditions ( p.l2| ) bound the yet unspecified integration constant r^, 
and, consequently, R^. 

Since ~ 10~^^, and assuming that /3, m ~ 0(1), one concludes that condition Eq. 
(3.12) is valid within the domain 1 kpc ^ ^ 50 kpc (where the flattening of the rotation 
curves is measured) for a wide range of values of r*: one could set, e.g., = 1 kpc. 

Regarding the validity of Eq. ( 3.13 ), the value r^, = 1 kpc is also sufficient to satisfy it 
within the domain 1 kpc < r < 50 kpc, since one expects tq ~ ^ 1 kpc and the (r/a)"^ 
factor (with a typically of the order of 1 — 10 kpc) is suitably suppressed by the v"^ = 10"^^ 
term. 

Hence, the assumption of a perturbative regime for the equations of motion Eq. ( ^.6| ) 
is self-consistent — although the deviation from geodesical motion is purportedly relevant, 
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dominating the Newtonian acceleration and leading to the flattening of the galaxy rotation 
curve. 

Two remarks are in order: firstly, one should note that the obtained expression for the 
non- minimal coupling f2{R) in the flattened rotation curve case /3 = is strikingly similar to 
that of Ref. where the non-trivial curvature term reads fi{R) = R[l + v'^ log(i?/-R*)] ^ 
R{R/R^,)'"°° . Following those authors, one may suggestively rewrite the geometric coupling 
with matter depicted in Eq. ( p. 101 ), so that action Eq. ( |2.1D reads 

'gd'^x, (3.14) 

where one defines a = v^/m or, conversely, one obtains the asymptotic velocity of the 
flattened region of the rotation curve through = ^Jam. 

Secondly, it is clear that this approach also displays the caveat of Ref. Q : assuming 
the same outer slope m for the visible matter density profile p of galaxies, the obtained 
flattening for the galaxy rotation curve is universal, that is, Uoo and /3 should be fixed pa- 
rameters of the model, unadjustable for different galaxies. This issue might be alleviated in 
a more realistic approach where one considers different outer slopes m for distinct galaxies, 
but it would be troublesome and somewhat artificial to reconcile models for the visible 
matter density with the corresponding asymptotic velocities and different outer slope 
behaviour of available rotation curves (that is, the variety of observed values for /3). 

One could conceivably invoke more complicated, yet unknown hidden dynamics within 
the model so to account for the variation from galaxy to galaxy — including different den- 
sity profiles, spatial symmetry and morphology. Nevertheless, this remains a noticeable 
disadvantage of this approach — and serves as an encouragement for considering alterna- 
tives. 

3.1 Logarithmically decaying galaxy rotation curves 

In the previous paragraphs, the issue of explaining the flattening of galaxy rotation curves 
by resorting to the non-geodetic motion arising from a non-minimal coupling between 
matter and geometry was addressed. For generality, this was performed assuming that the 
reported flattening does not imply on an actual asymptotic behaviour of f (r), and that 
inverse power law decays f^(r) oc are allowed. In this subsection, one further considers 
the possibility of a logarithmic decay in the outer region of the galaxy rotation curves, of 
the form of 




(3.15) 



As can be easily seen, this leads to an enclosed "dark matter" distribution of the form 
Mdm{r) oc u^(r) x r oc log(r/r*), so that the corresponding density is given by 

^TTr^Pdmir) = M^„(r) ^ paUr) « r'^ (3.16) 



S 



SttGR + ( ^ 
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This scenario is significant, as the obtained dark matter density coincides with the outer 
profile corresponding to the NFW model, which shall be object of further scrutiny in the 
following sections (see also Appendix A). 

Following the steps outlined in the previous paragraphs, one writes 



A/^(p(r)) = ^^logf-) 



(3.17) 



with solution 



A/2(/5(r-)) = -v^ 



2 ' * 



1 + loe 



As before, one sets A 



10 



-12 



Substituting the expression r 



(3.18) 

a a/poTp and the 



classical approximation R ~ SirGp into Eq. (3.18|) yields 



l/m 



1 + — log ^ 
m \ n 



(3.19) 



The perturbative conditions Eqs. (0) thus read 



m 



3 + log ( ^ 



1 + log 



ro\' 



< 1, 



(3.20) 



Following the discussion of Eq. ( p.l2| ) and (3.13), it is straightforward to check that these 
conditions are satisfied within the domain 1 kpc < r < 10 kpc for a wide value of r*, i.e. 
r* = 1 kpc. 



4. Pressureless dust with non-minimal gravitational coupling 

In the following, one assumes that the motion is fairly close to geodesical, posit an ap- 
propriate non-minimal coupling /2 (R) and ascertain under which conditions a dark matter 
component can be mimicked. As in the previous session, one isolates the effect of a non- 
minimal geometric coupling with matter by considering the usual scenario fi{R) = 2kR. 
Furthermore, one assumes that the effect of the non-(covariant) conservation law Eq. ( p.5| ) 
is negligible, so that test particles follow geodesies described by the metric gi^^y. This 
possibility will be scrutinized later on. 

One ultimately aims to establish a relationship between the behaviour due to dark 
matter and the one arising from the model arising from Eq. (^]^). Moreover, one aims 
to do so by mimicking accepted density profiles for the dark component of a galaxy. This 
reflects the overwhelming evidence towards some yet unexplained phenomena that behaves 
as non- or weakly interacting matter (asides from gravitational), and the willingness of 
the authors to account for those results, instead of simply conjuring another hypothetical 
candidate among the so-called "dark matter zoo" of possibilities. This marks a departure 
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from previous efforts, which also imply the flattening of the rotation curves, but cannot be 
interpreted in terms of an actual dark matter density profile. 

For generality, one considers a power law gravitational coupling of the form 



A(ii) = (!)". (4.1) 

with n a yet unspecified exponent and Rq a characteristic scale of the non-minimal coupling, 
given this choice for /2, one freely sets A = 1, so that there are only two model parameters, 
n and Rq. 

Notice that Eq. ( [4.lD does not exclude a more elaborate form for the non-minimal 
coupling: if one assumes that it is written as a series 

summed over the rational exponents p of the "full" model (integer or not), then the choice 
of Eq. (|4.l| ) may be regarded as an approximation valid where the p = n term is dominant. 
Likewise, this study is not in opposition to the findings of Ref. [^], which approached the 
possibility of including a linear coupling f2 = R/Ri — and, through the calculation of its 
impact on solar observables, it was found that Ri » SirGpQ, with pq the central density 
of the Sun. 

One considers that matter within a galaxy is modeled by a pressureless perfect fluid, 
i.e., dust; this is characterized by a density p, so that the energy-momentum tensor becomes 

Tf,u = pU^U, ^T = -p. (4.3) 
Using the Lagrangian density Cm = —p [14], Eq. (2.6) reads 



n 



^Ra) R 



R 



n 



R_ 

To 



p_ 

R 



+- 



2k 



and Eq. (p. 7]) becomes 



R 



1 

2^ 



1 + (1 - 2n) 



R 
Ro 



3n 



P 



-U 



i- 

\Ro 



p_ 
R 



(4.4) 
(4.5) 



The added terms induced by the non-minimal gravitational coupling f2{R) may be 
interpreted as contributing to the dark matter halo density profile poM- To see how this 
mimicking behaviour may arise, one first assumes that the non-minimal terms dominate 
very far from the galactic core (where the rotation curve flattens). Hence, in this exterior 
region, one may rewrite Eq. ([4.5D as 



4.1 Static solution 



Since the visible matter density p is assumed to be known, solving the above equation 
amounts to obtaining the relation R = R{p). It is easy to verify that the classical iden- 
tification R = 2kp fails, since one would obtain a differential equation for p, which might 
not be satisfied by the considered visible matter density profile. By the same token, one 
concludes that a solution defined implicitly by the expression 



l-2nfRY 



(4.7) 



is self-consistent, since the gradient term vanishes everywhere; for this reason, one dubs 
this as a "static solution" . This yields 



R — R\ 



l-2n p 
2k R^ 



1 1/(1-") 



R(] 



(1 



2n)^ 
Po. 



1 



(4. 



introducing the characteristic density po = Rq/2k,. 

Recall that the objective of this study is to obtain an effect that mimics a dark matter 
profile pdm- One is tempted to extract this directly through pdm = 2k,R; this implies 
that the dominance of the non-minimal gravitational coupling f2{R) effects is equivalent 
to the condition pdm > P, valid when the galactic rotation curve flattens. However, this 
identification assumes that the effect of the non-minimal coupling f2{R) may be interpreted 
as a pressureless dust distribution; to qualify this, one should first analyze the full field 



Eqs. (4.4), assuming solution Eq. (4.8) and a perfect fluid distribution for the dark matter 
component one wishes to identify, namely 



T{dm)fjLV = (Pdm + Pdm) VfiVp + PdrnQfiu, 



(4.9) 



where pdm is the density, pdm the pressure and the four-velocity of the mimicked dark 
matter component. 



R 



1 - 2n 



R 



I - 4n 

1 R \ 



(4.10) 



[{Pdm +Pdm) Ufj^Uu +Pdm9tiu] ■ 



By assuming that the four-velocity of the mimicked dark matter component is equal to 
that of ordinary matter, = U^, one obtains the pressure and density of the former: 



Pdn 



and, substituting Eq. (W^, 



1 — n 
1- in' 



-2kR 



n 



Pdn 



1 - 4n 



2kR, 



(4.11) 



Pdn 



1 — n 
1 - An 



Po 



[1 - 2n) 



P_ 



1/(1-") 



(4.12) 
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Thus, one concludes that, sufficiently far away from the galactic core, where the rota- 
tion curve has flattened, the effect of the non-minimal gravitational coupling with ordinary 
matter mimics a dark matter component which is dragged by ordinary matter and amounts 
for a non-vanishing pressure and an equation of state of the form 

Pdm = Z Pdm- (4.13) 

1 — n 

This dark matter-like component may have a negative density pdm < 0, if 1/4 < n < 1; 
although strange, this is not a pathology of the model, as the curvature in the considered 
outer region will remain positive, i? > 0. 

The above treatment was pursued with no regard concerning the problematic values 
that the exponent n might take on. Indeed, one sees that n = 1 produces a matter 
distribution with no density and negative pressure, while n = 2 leads to a zero curvature 
solution i? = (see Eq. (|4.8|) ) and n = 1/4 yields a sing ular Eq. (p^) . A correct 
treatment of this cases is not developed here, as one shows in the following section that 
adequate dark matter density profiles may be mimicked by the non-minimal coupling with 
an exponent n different from the discussed values. 

4.2 Mimicking dark matter density profiles 

In the context of dark matter models, it is assumed that, at the region where the galaxy 
rotation curve has flattened, both visible and dark matter might be modelled by a power 
law density profile, characterized by the outer slopes m and m', respectively (see Appendix 
A for a thorough discussion): 

Pv iy-j , Pdm ^ Pdi— j , (4.14) 

so that, in this outer region, they are related by a power law of the form pdm oc pP, defining 
the slope ratio p = m! jm. The constant /O^, is specific to each visible matter distribution 
p, and will be related to this quantity and to the previously defined /jq- 

This relation is crucial to our study, since a similar power law was obtained in the 
previous session: one now resorts to Eq. ( [4.12 ) to provide a translation between the 



parameters -Rq j ^ of the non- minimal coupling model and the mimicked dark matter density 
profile quantities a', m' and p^, assuming a given visible matter density profile characterized 
hy pv, m and a. 

Firstly, one obtains the following relationship between the exponent n and the outer 
slopes m and m', 

1 m' ml , , 

= — = p->n = l = 1 , (4.15) 

1 — n m m p 

which yields the power-law discussed above, 
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implying that 



Since, in this study, one assumes that dark matter is a dynamically generated effect arising 
from the non-minimal gravitational coupling /2 (R) , the length scale a' is inherited from the 
visible matter density, so that a' = a. However, one should note that the inner behaviour 
of our model no longer has Eq. (^^) as a solution, since one can no longer rely on the 
condition p^m > P- This shall be discussed in the subsequent session. 

4.3 Discussion 

In order to verify the consistency of the mechanism proposed in the previous section, one 



now ascertains the behaviour of the trace Eq. [4.5| , without the assumption of solution Eq. 
( |4.8| ) nor the outer region approximation pdm ^ P- To do so, one first rescales all relevant 
quantities; one considers the length scale rg = and the already used characteristic 

density po = 2k,Rq (so that, if n < 0, GR corresponds to Rq = 0, that is, tq — )■ oo). These 
define the dimensionless quantities 



a = - r = - = ( l.Y^^^"'^ g^^-lJL (418) 

ro ' '^o ' \PoJ ' ~ PqO Ro' 

and one also defines the dimensionless visible matter density 0* = p/ Pq = for conve- 

nience. Introducing the above into Eq. (4^), one obtains 

g = e-^ + {l-2n)d' -^^Ud"\ (4.19) 

a 

where □ = TqD = \I\/Ro is the rescaled D'Alembertian operator. In this form, the classical 
identification R = —2Kp reads q = 0~^, while the solution Eq. (|4.8| ) becomes q = {\ — 

2„)l/(l-n), 

For simplicity, one assumes that the galaxy is Newtonian, that is, that one may disre- 
gard relativistic contributions to the gradient term in the above equation. This said, one 
writes the D'Alembertian in spherical coordinates, assuming the spherical symmetry im- 
posed by the visible matter density p{r) (it is a simple matter to extend this to cylindrical 
coordinates and axial symmetry): 

dp 2 d ,^ , 

□ =-2+3-. (4.20) 

ar^ r ar 

Introducing the inverse coordinate y = = rg/r, one obtains the simplified form 



□ = 2/'^' (4-21) 

so that Eq. ( p^ becomes 
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Q = e-^ + {l- 2n)Q^ - ^y^-^. (4.22) 



6n 4dV"^ 

One cannot straightforwardly extract the GR result q = from the above equation, as 
there is no model parameter allowing one to take the appropriate limit. Notice, however, 
that in GR Rq = Q and ro diverges, so that the equation vanishes trivially. 

The programme consists in introducing the visible matter density 0*(y) = 0^~"'{y) into 
the above, and extract the rescaled curvature g, from which one might infer the effect of 
the non- minimal coupling f2{R), as well as the validity of solution Eq. (|4.S| ). To do so, 
one resorts to the Hernquist density profile (with inner slope 7 = 1 and outer slope m = 4, 
see Appendix A), written as 



Pir) = --——3' (4-23) 



Ma 1 
2vrr {r + af 

so that = M/{2'Ka^) in Eq. ( [4.14 ), where M is the total visible mass of the galaxy. This 
yields 



„,4/(l-n) „,m' 

e{y) = (2(5a)i/(i-")— 7^ = {25(1^''^- ^ ,,, , 

^ ' (1 + ay)3/(i-'^) ^ ' (1 + ay)3'"'/4 ' 

using relation ?i = l — l/p=l — 4/m' and defining the dimensionless parameter 

i = ^ = !i, (4.25) 

with Tg = 2GM/(P' being the Schwarzschild radius of the galaxy. Clearly, in the large r 
limit 6 oc y™" is proportional to the desired dark matter density profile. 

In order for the assumption of negligible general relativistic corrections to be self- 
consistent (expressed in the flat spacetime form of the D'Alembertian operator, Eq. ( [4.20| )), 
one assumes that the visible matter density is sufficiently low, so that these corrections do 
not dominate the effect arising from the non-trivial gravitational coupling. This translates 
into the condition (5 < 1, which shall be tested later. 

4.3.1 Characteristic density and background matching 

One can now infer the characteristic density p^: 

where poo = SHq/^SttG) is the critical density of the Universe, and Hq ~ 70 {km/ s) /Mpc 
is the Hubble constant. Notice that, in a cosmological setting, the non-trivial coupling 
(R/Rq)"' is given by (poo/po)"', so that the obtained relation Eq. ([4.26| ) might support a 
previous claim: that the proposed form f2{R) = (Rq/R)^ could hint at a unified description 
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of dark matter (as presented here) and dark energy, responsible for the acceleration of the 
expansion rate of the Universe, which is relevant in a cosmological context. 

This avenue of research is not followed here, but one should nevertheless point out 
that this might account for a feature of the mechanism described in this work: not only it 
gives rise to a "dark matter' component with adequate density profile, but it also endows 
it with negative pressure Eq. ( [4.13| ) (for n < 0) — an attribute usually associated with 
dark energy. 

This feature is not dependent on the particular geometry here adopted, and is not 
exclusive to the proposed power-law coupling f2{R) (with negative exponent). Indeed, 
it can be shown that the "curvature pressure" associated with the effect of either fi{R) 
or /2(-R) involves terms in the corresponding derivatives Fi{R) and F2{R) (see e.g. Ref. 
fl^ for a proposal of usual f{R) theory in a cosmological context). As a result, the non- 
minimal coupling scenario yields a negative "pressure" if F2{R) < 0: in particular, if it 
takes a power-law form f2{R) oc R" with a negative exponent n. 

Hence, as the visible matter density p asymptotically falls to its background cosmo- 
logical value, so does the dark matter component pdm, thus presenting a possible smooth 
transition from the local galactic scenario to the averaged, large scale cosmological setting. 
Other models of unification of dark energy and dark matter include the Chaplygin gas 
model [17, |l^, |2^] and some particular field theory constructions pl|]. 

This may be used to signal the distance r^o where the total galactic density profile 
(including dark matter) blends with the averaged cosmological value, defined by 2KR{roo) = 
Poo- Assuming the "static" solution g = {1 — 2n)^/^^~'^\ this yields 



%oo) = ^(l-2n)-Va-), 
Po 



which, in the large r limit, reads 



(4.27) 



(2<5a)i/(i-")yi/(i-") = ^(l-2n) 

Po 



-l/(l-n) 



2(1 - 2n)6a 



3 X 3"(1 



2n) 



1/4 



roH 



(l-n)/2 



r.r, 



1/4 



(4.28) 



One can interpret the distance Too as an indicator of the radius of the mimicked "dark" 
matter halo; for consistency, it should vanish when the non-minimal coupling Eq. (J^^ is 
negligible, i.e. i?o — ^ if the exponent n is negative or i?o ~^ for n > 0; this is trivially 



satisfied, since Too oc r^^'^ = i2g " — (excluding the constant coupling n = 0). 



4.4 Crossover between visible and dark matter dominance 



The key feature of the present proposal is that the behaviour of the non-minimal gravita- 
tional coupling, which enables the scaling law Eq. ( [4. 12]) , (cf. below Eq. ( 4.18| )) yields an 
appropriate dark matter-like component from a given visible matter density profile. The 
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model parameter Rq introduces both the length scale rg as well as a density scale po; given 
that the discussion of dark matter is related to the flattening of the galaxy rotation curve 
in the outer region, where it is assumed that ro has a clearer physical meaning — it should 
signal the crossover distance Tc between the classical solution R = 2kp q = and the 
outer solution R = Rq [(1 - 2n)(p//>o)]^''^^~"^ ^ g = {1 - 2n)i/(^-"). This should occur 
when both solutions are similar, 

0-" = (1 - 2n)^/(^-") ^ e^y) = (1 - 2n)-i/". (4.29) 

Assuming that this crossover occurs about the outer region of the galaxy, so that Tc > a, 
one might still approximate the visible matter density profile by its outer slope profile, 
obtaining 



e^y) = 25rl-^ = (1 - 2n)-'/^ ^ (4.30) 
= 2(1 - 2n)^/"5r[ja = 2(1 - 2n)i/"r,rga. 

Hence, the crossover radius Tc is not simply proportional to rg, as one might naively 
expect, but involves a sort of "geometric mean" of all the length scales present. This 
translates into the dimensionless crossover coordinate yc, given by 



2(1 - 2n)^/"5a. 



(4.31) 



4.5 Behaviour of the gradient term 
4.5.1 Inner region (r < a < r^.) behaviour 

The discussion above assumes that the gradient term on the r.h.s of Eq. ( [4.22| ) is negligible 
— which is valid for the dark matter dominated region, g ~ (1 -2n)i/(i-"). However, one 
should analyze the relevance of the gradient term present in Eq. ( [4. 22]) as one approaches 
the inner regions of the galaxy, y — )• oo, since in this region the power-law Eq. ( [4.14 ) is 
no longer an accurate description of the visible matter density profile, and large gradients 



may arise. From Eq. (4.22), one sees that the gradient term is negligible if 



6n- 



e dy^ 



6n 



y4 d^9^ 



dy^ 



At distances r < a, the Hernquist density profile reads 



< 1. 



(4.32) 



p{r) 



M 



0*{y) = 25- 



;2 ■ 



(4.33) 



271 a? r 

In the interior region where the visible matter density is ruled by its inner slope be- 
haviour, Tc < a ( dye > 1), the consistency condition Eq. ( 4.32| ) reads 



6n- 



7^ d^^r 



dy' 



Qn\n - l)y"+^ ( ^ 



26 



<. 1. 



(4.34) 
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One first studies the n > — 1 case, thus defining the quantities yk{in) ^-nd r^j-j^) 



y < Vkiin) 



1 



6n^(l — n) \a 



25 



l/{n+l) 



(4.35) 



6n^(l — n) 



2 \ 1^" 

2r. 



l/(n+l) 



As one approaches the center of the galaxy, y — )• oo, the relevance of the gradient term of 



Eq. ( 4.22 ) increases, so that eventually one gets y > yfc(j„)(at shorter distances one would 
have to include relativistic corrections of the form Vg/r, which are not considered here). 

Although not exhibited here, numerical results show that this behaviour is not a direct 
consequence of the visible matter inner density profile: instead of the inner slope of the 
Hernquist distribution, one might choose a density profile that flattens towards the galactic 
core (that is, so that {dp/dr)\r=o = — see e.g. Ref. ^^); however, since \dp/dr\ 
increases outwards, the gradient term in the r.h.s of the equation of motion Eq. ( |4.22| ) 
would eventually come to dominate, albeit at a different region. The bottom line is that, 
for any given visible matter density profile p, the gradient term will eventually dominate 
below a threshold r < r;j(j„) . 

However, the argument implicitly assumes that the Newtonian regime extends down to 
these core regions; although one does not pursue a full treatment of the relativistic regime 
(when r < r^), one expects the behaviour of the gradient term to be significantly different 
in these region, since the D'Alembertian operator is no longer of the form y^d^/dy^ (in 
particular, corrections of the form r^/r = 6y arise, which are deferred to a future study). 
Hence, one may circumvent this issue by ascertaining that rfc(j„) < — t- > 1. This 

yields the condition 



;1— n 



,1— n 



In the case n = — 1, condition Eq. ([4.34|) reads 



Tgro 



> \/3, 



(4.36) 



(4.37) 



which coincides with Eq. ( [4.36 ). 



4.5.2 Tracking behaviour of gradient solution 

As discussed above, a suitable choice of the model parameter ro ensures that the gradient 
term on the r.h.s. of Eq. ( [4.22 ) dominates only within the relativistic region r < rg, where 
the Newtonian approximation breaks down. However, a numerical analysis (which will be 
displayed later, for two relevant cases) indicates that it is impossible to further tweak the 
model so that this gradient term will not become dominant after the cross-slope region 
r ~ a, while still mimicking the dark matter behaviour. Hence, one concludes that the 
"static" solution Eq. (|4.8[), given by ^ = (1 — 2n)^/(^~") is not valid: in particular, the 
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discussion of the dominance of the gradient term in the outer region a < r < rc is flawed, 
as the crossover distance Tc assumed for its vahdity. 

Hence, instead of the crossover Vc, one must consider the distance rk> a above which 
the gradient term dominates (or, equivalently, a dimensionless coordinate yk) < one 
may write it as = — )• = l/(Ca) (analogously to the relation rc = ea used to study 
the "static" solution). 

Using the outer slope Hernquist profile, 0*(y) = 25ay^, Eq. ( 4.32| ) reads (for negative 

n) 



6n 



?/4 (fO'^ 



24n2(l - 4n) (25a)""^ y^"-^ ^ (4 gg^ 



so that, when y = yk = 1/(C'^)) one has 

12 X 2"n2(l - 4n)5"~iai-3"^2-4n ^ (4 39^ 

which is used to set tq, for a given n. 

One must now look at the behaviour of Eq. ( 4.22| ) when it is dominated by its gradient 
term, thus reading 

fl^^ 0_^^^_^Say^;f^^^ 

dy"^ 6ny^ 6n ' 

having substituted the outer slope Hernquist profile, 0*(y) = 25dy'^. This equation admits 
what one dubs as a gradient solution: clearly, this cannot be a constant solution q ~ const, 
(like the "static" one Q = {I - 2n)^/(i-")). 

Since one is interested in the behaviour of the above equation in the outer region r > a, 
it is tempting to resort to a series expansion around y = 0. However, this would require 
the stronger assumption r > tq, which may not be valid; in particular, if the characteristic 



density po has an order of magnitude similar to the cosmological background pQ, Eq. (|4.2q ) 
will yield tq ~ 100 Gpc, which is much larger than the length scale of a galaxy. However, one 
may instead define the coordinate z = ay = a/r and the dimensionless function ijj = q^~^, 
rewriting Eq. ( 4.40| ) in the Emden- Fowler form 



^ = 0^4„/(l-„)^l/(n-l)^ (4.41) 



with 



G = — i . (4.42) 

For a given n, the solution to this differential equation may be obtained implicitly, in a 
parameterized form z = z{t, G, n), = g{T, G, n)) instead of pursuing this possibility 
or analyzing some relevant cases (such as the n = —1 and n = —1/3 scenarios studied 
in a following section), one undertakes a more physically oriented, albeit less rigorous, 
exploration of the values of interest. 
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In order to do so, one may ascertain the order of magnitude of the r.h.s. of Eq. ( 4.4C| ) 



as the gradient term begins to dominate, at y = yk = £,a. One considers the classical GR 
solution g = valid up to this point, that is 

d'g''-'\ _ e'-{yk) _ eM _ Sa ^^^^^ 



dy'^ Jy=y, ^nyl 6nyl 3n' 

For typical values of a = 1 — 10 kpc and = 10~^ — 10^^ pc (corresponding to galaxies 
of approximately lO^'' — 10^^ solar masses), condition Eq. ( |4.36D yields > 1 Gpc, much 
larger than and a; hence, one concludes that the r.h.s. of the above equation is much 
smaller than unity, independently of ^. After the gradient term of Eq. ( [4.22 ) begins to 



dominate, the dimensionless quantity g is essentially ruled by the differential equation 

and, since g decreases, one concludes that the gradient solution to Eq. ( [4. 221 ) is approxi- 
mately constant when ay — )• — and is similar to the "static" solution ^ = (1 — 2n)^^^^~'^\ 

4.6 Tully-Fisher law 

The Tully-Fisher is an empirical relation between the intrinsic luminosity L and a power 
of the rotation velocity Vqo of a spiral galaxy (with the related Faber-Jackson relation for 
elliptic ones). Since the luminosity is proportional to the visible mass, L oc M, it may be 
written as M oc v'^, with 3 < o" < 4 for spiral galaxies. 

The velocity width Voo, measured in the outer region of a galaxy, depends upon the 
total mass distribution within the delimited zone. This may be obtained by integrating 
the profile 2kR = PqOq, after numerically solving Eq. ( [1.22| ); recall that, from Eq. ( 4.11| ), 



this also includes the effect of the negative pressure arising from the non-minimal coupling 



/2(-R). Clearly, Eq. ( 4.12 ) indicates that the mass profile of the visible and dark matter 
components will differ substantially. Actually, the total mass corresponding to the NFW 
dark matter density profile is infinite. 

This said, one can still get a crude estimate of the form of the power-law relating M and 
Voo, deferring a more elaborate numerical estimate to the subsequent session. One assumes 
that M(im{f) oc r at large distances (leading to a flat rotation curve), and Mdm ^ M. This 
yields the familiar oc Mdm relation that, resorting to the scaling Eq. ( [4.12| ), yields 

Mocw2(i-"). (4.45) 

For the scenarios considered in the following session, one obtains M oc v^^^ (for n = 
-1/3) and M oc (for n = -1). 

4.7 Conservation of energy-momentum tensor 

In section ^, one has first attempted to justify the flattening of the galactic rotation curves 
as due to the non-conservation of the energy- momentum tensor of matter, i.e., through 
the non-geodesical motion. However, one is now exploring a mimicking behaviour of the 
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model under scrutiny, so that a dark matter-like behaviour emerges as a consequence of the 
dynamics of Eq. ( 4.22| ). Hence, one is assuming that matter follows geodesies, as explicitly 
stated in the beginning of section After the previous characterization of the model, one 
is able to ascertain the validity of this assumption. 

One begins by rewriting Eq. ( |3.2D , with p = for pressureless matter, the Lagrangian 
density is taken to be £ = —p and the adopted non-minimal coupling Xf2{R) = (R/Rq)'^: 



with radial component 



-n- 



1 + 



r 



-n- 



log( 



log 



R{r) 

Rn 



(4.46) 



(4.47) 



assuming spherical symmetry and the Newtonian approximation hrr- = grr- — {vr/c)'^ ~ 1- 



4.7.1 Non-physical result arising from the "static" solution 



From the discussion presented in subsection 4.5, one recalls that the gradient term in Eq. 



(4.22) will become dominant at distances r > r^, so that the crossover between visible and 
dark matter dominance does not occur at the distance Tc discussed in subsection [4.4] . This 
is not only a numerical mishap, but a crucial component of the model presented here: if 
the "static" solution g = (1 — 2ny^^^~'^^ is valid for the outer region r > a, the dominance 
of dark matter would be equivalent to the condition (1 — 2n)f)" ^ or, from definitions 
Eq. ( p^ , {R/RoT > 1/(1 - 2n): if this was the case, Eq. (|446|) would read 



r 



-n 



-n- 



R!_ 



(4.48) 



Recall that, in the outer region r > a where the gradient solution to Eq. (4.22) tracks 
the constant "static" solution g = (1 — 2n)^/^^~'^\ so that R oc r~"^' follows the desired 
dark matter density profile. Hence, the above equation would lead to a force fr = nm' /r 
(pointing towards the center of the galaxy, since n < 0). Clearly, for small \n\ and m', this 
would lead to a flattening of the rotation curve on the onset of the outer region r > a, 
with a asymptotic velocity Voo/c = \/ —nml ~ 1, in blatant violation of the observed value 
Woo ~ 100 km/s ~ 10~^c; asides from this non-physical result, it would undermine the 
primordial objective of the proposed programme: to obtain a dynamical effect from the 
non-minimal coupling that mimics dark matter, so that geodesical motion produces the 
flattening of the galaxy rotation curve. 



4.7.2 Geodesical motion from the "dynamical" solution 

As discussed above, the crossover between visible and dark matter occurs not due to the 
onset of condition {R/Rq)^ > 1, but because the gradient term in Eq. ( [4.22| ) becomes 
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dominant. Hence, one can instead take (R/Rq)" < 1 and consider for Eq. ( 4.47] ) the 
approximation 



r 



-n 



RoJ 



n-l 



Rq 



(4.49) 



Clearly, the pathological result stemming from Eq. (4.4S) is suppressed in this case. Aiming 
at comparison with the already defined dimensionless quantities, one rewrites it as 



r 



-n{Qer-\eO)' 



An 



1 



n 



{25a) 



i/(i-n) <f_ (rQ_Y^/{^-n) 



c 
r 



(4.50) 



after using 6 = 6\^^^ "\ the outer Hernquist density profile 9^ k, 2<5a(ro/r)^ and the asymp- 
totic result Q ~ const.. In order for this force to dominate the Newtonian gravitational 
force /at = —GMjr^ = —{rs/2){c/r)'^ (which is smaller than the overall gravitational force, 
including the dynamically generated dark matter), one must have 



8n 
1 — n 

r > ro 



{25af 



/(l-n) |^royn/(l-n) ^ 

r 



1 



n 



8n 



r J 

l-n 



(4.51) 



a 



l/(l-5n) 



For the already considered values of a = 1 — 10 fcpc and = 10^^ — 10^^ /cpc, one 
concludes that it suffices that ro > 200 k-pc for the above lower bound to be larger than 
1 Mpc, thus ensuring that the extra force Eq. ( [4.49 ) arising from the non-conservation of 
the energy-momentum tensor is negligible at galactic scale. 

4.8 Energy conditions 



Following Ref. [^^, one discusses the relations required to satisfy the strong, null, weak, 
and dominant energy conditions (SEC, NEC, WEC and DEC, respectively); one considers 
p = 0, /i = 2kR and /2 = (R/Rq)'^; using the equality 



f2 



R 



2np 



1 + 



(4.52) 



the energy conditions read j24] 



1 + 



P>0, 



(4.53) 



with 



E„ 



■ -2n 


, SEC 





, NEC 


< 

2n 


, DEC 


. n 


, WEC 



(4.54) 
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Clearly, both the NEC and SEC are always satisfied, for a negative power law n < 0. 
The WEC and DEC are also fulfilled for the ranges — 1 < n < and —1/2 < n < 0, 
respectively. 

Moreover, as will be seen, the condition R > Rq is always valid at galactic scales (due 
to the dominance of the gradient term of Eq. ([4.6|)) for the considered models n = —1 
(isothermal sphere profile) and n = —1/3 (NFW profile), so that (1 + {R/Rq)~"')~ 0; 
hence, the WEC and DEC are also satisfied for the isothermal sphere and NFW profiles. 

Asides from these energy conditions, one should also verify whether the considered 
model does not give rise to the so-called Dolgov-Kawasaki instabilities, whereas small cur- 
vature perturbations could arise and expand uncontrollably, thus rendering the underlying 



theory unphysical ||2^. From Refs. |24, 26 1, one writes the mass scale m^)/^ for these 
instabilities as 



2 _ 2k - im) + f2iR)R] + mm _ 



2kR(^) "'+n{2n-l)p 

—R, 



2n{n — l)p 

using T = Cm = —p- Since the non-minimal coupling mimics an additional, positive 
matter component, one has p < 2kR; furthermore, the cosmological matching condition 
Eq. ( [1.26| ) hints that the characteristic density obeys 2kRo = pq <^ p, so that one may 
write R/Rq ^ 1. Since n is a negative exponent, one may approximate the above equation 

by 

.-(^] ^ . (^] - ,4.50) 



-Ro / n(n — l)p \Ro J p 2n{n — 1 
which, in terms of the dimensionless quantities defined in Eq. ([4.18|), reads 



mlK^-Q^'^'O (4.57) 
2n{n — 1) 

Since Tn'j-,j^ is negative, there is a potential concern about instabilities. These evolve at a 
length scale rj^x 

^ J-^ro = V2n(n-l)£>"-i/2ri/Vo. (4.58) 

R / \ Po 

At large distances, when the dark matter-like component dominates and the scalar curva- 
ture is larger than the matter density, the dimensionless quantity g approaches a constant 
value (displaying the tracking behaviour of the "gradient solution", as discussed before). 
The remaining quantity 9 behaves proportionally to the mimicked dark matter distribution 



with an outer region profile (see Eq. ( 4.24 )), 



9ir) = {26ar'/\a/rr', (4.59) 
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so that one may write 



roK = VMn-l) ..Vw/« - ^0. (4.60) 



(2(5a)™'/8 

Thus, one finds that the typical length scale r^x grows linearly, for the m' = 2 isothermal 
sphere scenario, or with r^/^, for the m' = 3 NFW density profile. 

The perturbative expansion leading to the mass scale Eq. ( 4.55| ) assumes an initially 



small perturbation 5R to a constant background curvature R; hence, one must check if the 
obtained length scale r£)K is larger than the characteristic distance over which one cannot 
assume that the curvature is constant. The later is given by L = —R{r)/R'{r) = r/m'. 

Comparing the instability length scale r^K against the constant curvature domain L, 
one concludes that the validity of the perturbative expansion, embodied in the inequality 
r£)K < L, requires 

fi"-V2 /r\m72 r 

iV'"-' < (,-«/2(2ia)"'/8 ! 



Y^2n(n — l)m' 

Inserting the scaling relation m' = m/(l — n) and the outer slope m = 4 of the visible mass 
Hernquist density profile, one rewrites this as 



l + n - , 



«y <^^-"+^/^(25a)^-y^. (4.62) 

This relationship shows that the perturbative expansion is invalid, in the case of the n = — 1, 
isothermal sphere distribution, or limited to a distance r < (where the Newtonian regime 
breaks down), in the case of the n = —1/3, NFW profile. Hence, the instability detected 
via Eq. ( |4.57D is not physically relevant. 

5. Mimicking the dark matter density profile of spherical galaxies 

One now aims to further develop the mechanism described in the previous sections, by 
attempting to model the rotation curves of galaxies of type EO, which exhibit an approxi- 
mate spherical symmetry adequate for the results obtained in this work. Since its manifest 
purpose is to avoid the intricacies related to the accurate modelling of galaxies, one opts 
for using seven readily available galaxy rotation curves, as provided in Ref. [|7|,|8|,||,||, 



divided into visible and dark matter components of the following galaxies: NGC 2434, 
NGC 5846, NGC 6703, NGC 7145, NGC 7192, NGC 7507 and NGC 7626 galaxies. 

In the discussion presented here, particular attention is paid to the scenarios whereby 
the proposed mechanism mimics either the NFW or the isothermal sphere dark matter 
density profiles (given by the n = —1/3 and n = — 1 cases, respectively), due to its 
ubiquity in the literature. Moreover, the analytical approach taken assumed that the non- 
minimal coupling was given by a single power-law f2{R) = {Rq/R)^^ thus yielding a natural 
explanation for either of the two profiles. 
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Conversely, it is known that some galaxies are much better fitted by the NFW profile, 
while others conform to the isothermal sphere model: hence, arguing for a single power- 
law coupling between matter and geometry amounts to defending a particular dark matter 
density profile — a task clearly outside the scope of this work, and better left to the 
astronomy community. 

Furthermore, the present framework offers a convenient way through which both dark 
matter density profiles may be implemented simultaneously, with no a priori need for any 
unnatural selection. One is thus led to conclude that the assumed single power-law coupling 
should be generalized to a sum over exponents n, hinting perhaps at the Laurent series 
expansion (since n is not constrained to positive values) of a more evolved form for f2{R)- 

Given the motivation above, one undertakes the numerical task of modelling the rota- 
tion curves of the selected galaxies by combining both the NFW as well as the isothermal 
sphere dark matter density profiles, i.e., choosing the non-minimal coupling 




/2(^) = y^ + ^, (5.1) 
which, by taking the trace of the Einstein Eq. (4.4), leads to an equation analogous to Eq. 
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Notice that one can no longer resort to the identification of a "static" solution, as before: 
this is due to the non-linearity of the problem, so that adding the n = — 1 and n = 
—1/3 non-minimal couplings does not translate into a sum of the solution obtained before. 
Likewise, one can no longer resort to the straightforward identification with a "dark" 
matter energy- momentum tensor endowed with an equation of state such as Eq. ( 4.13| ). 
Tentatively, one can posit that the effects of combining the two coupling amount to adding 
the corresponding mimicked "dark" matter density profiles, together with interaction effects 
arising from the mentioned non-linearity. 



One now attempts to define the dimensionless quantities analog to Eq. ( 4.18 ), aiming 
at a generalization of the equation of motion Eq. ( [4.19| ); since there are two lengthscales 
ri = R^^ and r^ = R^^, one could adopt a normalization lengthscale ro = /(?"i,?"3), so 
that ro —7- ri when — t- +oo ( i.e., the R3 coupling is "switched off"), and vice-versa. 
However, this would unnecessarily complicate the obtained equation of motion, specially 
since one may isolate one of the power-laws in the non-minimal coupling Eq. (^.ip. 
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This said, one adopts the simpler normalization ro = ri, thus defining the dimensionless 
quantities as in Eq. ( 4.18| ). Some algebra yields the equation of motion 
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with the parameter e = {R^/ RiY^^ = (ri/rs)^/^ indicating the relative strength of the two 
couphngs present in Eq. ( |5.1| ): the pure n = —1/3 NFW scenario is obtained when e — )• cxd, 
while the n = — 1 isothermal sphere profile is given by the limit e — )• 0. 

The programme for the numerical session is simple: for a given rotation curve, one 
uses the visible matter density profile p as an input to Eq. (^^), adjusting the parameters 
ri and so that the obtained full (visible + dark) matter density profile yields the best 
fit available. Ideally, the model parameters ri and r^, are universal, so that all rotation 
curves should be fitted simultaneously. However, one instead proceeds by obtaining pairs 
for these quantities corresponding to each galaxy, and then discusses the obtained results: 
this allows for the identification of any abnormal case, possible trends, etc.. 

Following the previous sections, one attempts to model the visible matter density with 
a Hernquist distribution. As it turns out, the visible components of the selected rotation 
curves are not compatible with just one of these curves, but are very well fitted by two 
Hernquist profiles: each is characterized by a mass Mi and lengthscale Oj (i=l, 2), with 
Ml < M2 and oi < 02 — thus modelling a core + diffuse matter distribution, with the 
latter dominating. 

The results obtained from the numerical solution of Eq. ( |5.4D with the described 
approach are presented in Table |l], which displays the values for Mi, M2, oi and 02; the 
best fit values for the characteristic lengthscales ri and are depicted in Table |2|. These 
are ihustrated by Figs.|l| and | (left): as can be seen, the composite non-minimal coupling 
Eq. (|5.1| ) provides close fits to all rotation curves in the outer region, with some galaxies 
exhibiting a discrepancy in the inner galactic region: this could be due to uncertainty in 
the original derivation of the rotation curves or enhanced effects due to the deviation from 
purely spherical symmetry (which should be more pronounced near the core). 

The density profiles are also plotted for visible matter and full visible + "dark" matter, 
together with the contribution arising from the gradient term present in Eq. ( |5.4| ) in Figs. ^ 
and H (right): the contribution from the latter is completely overlapped with the mimicked 



"dark matter" profile, indicating that it dominates the r.h.s of Eq. (5^). Although not 
presented, for brevity, it is straightforward to check that this dominance is supported by 
the criteria discussed in section [4.5| ; this verification is achieved by resorting to the naive 
identification a = oi, M = Mi ( i.e. considering only the long-range Hernquist profile that 
dominates the distribution of visible matter) — and disregarding the non-linear effects 
arising from the interaction between the two power-laws present in the composite non- 
minimal coupling Eq. ( ^.1[ ). This extrapolation attempts to reconcile the analytic results 
obtained for the simpler scenario addressed before the numerical fitting, and constitutes a 
qualitatively cross-validation of both efforts. 



5.1 Single power-law non-minimal coupling 
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Table 1: Relevant quantities for the selected set of galaxies: AI is the visible mass, Ts the cor- 
responding Schwarzschild radius, ai, 02 and Mi, M2 the best fit values of the two-component 
Hernquist profile of the visible matter density. 
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Table 2: Best fit values for the characteristic length- 
scales ri and ra for the composite and separate fits of 
the galaxy rotation curves with the n — —1 isothermal 
sphere and n = —1/3 NFW mimicked "dark matter" 
scenarios, together with background matching distances 
roo 1 {n — —1) and r^o 3 {n — —1/3) for the composite 
non-minimal coupling, 00 in the composite scenario 
indicates that the corresponding scale Ri = . Units: 
ri (Gpc), rs (10^ Gpc), Too 1 and r^o 3 (kpc). 



In this paragraph the results of 
fitting a single power-law non-minimal 
coupling f2{R) = (i?/i?o)" (with n = 
— 1/3 or n = —1) to the rotation 
curves of the selected galaxies are pre- 
sented. Table ^ indicates the val- 
ues obtained from the best fit scenar- 
ios for the characteristic lengthscales 
ri and r^; as in the composite cou- 
pling Eq. {5A), these are of the order 
ri 10 Gpc and ~ 10^ Gpc. The 
obtained fits, although of worse qual- 
ity, exhibit a similar variation: an 
average fi = 16.8 Gpc with a stan- 
dard deviation ai = 10.5 Gpc and 
fs = 1.45 X 10^ Gpc and 0-3 = 1.26 x 
10^ Gpc. 



The corresponding rotation curves 
are depicted in Figs. ^ and |5|; as can 
be seen, the composition of the n = — 1 and n = —1/3 scenarios (Figs. |l| and |2|) improves 
the quality of fit considerably. For illustration purposes, the visible matter and full density 
profiles are depicted for the NGC 7145 galaxy in Fig. ^; the density distributions arising 
from the gradient and static solutions are also shown in these figures, with the tracking 
behaviour of the latter visible beyond the crossover between visible and dark matter dom- 
inance, for r > Tc ~ 10 kpc. This galaxy is chosen because it offers fits of similar quality 
with both the NFW and the isothermal sphere models. 
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5.2 Characteristic density and background matching and mimicked "dark mat- 
ter" inter-dominance 

As can be seen, the rotation curves of the selected galaxies are fitted by values for the 
characteristic lengthscales of the order ri ~ 10 Gpc and ~ 10^ — 10^ Gpc. Keeping in 
mind that the analytic results obtained in the previous sections concern a single power-law 
non-minimal coupling only, not the composite form Eq. ( |5.1| ), one may nevertheless use 
Eq. ( [4.28| ) to obtain a quantitative figure for the background matching distances. 



^x3"(l-2n) 



1 1/4 / X (l-n)/2 



3.8 Gpc \ 1^ ^2„^l/4 



r.^rfa) , n = — 1 



after inserting the Hubble constant Hq = 70.8 {km/ s) /Mpc. Similarly to the caveat above, 
one also notices that this expression was obtained with the assumption of a single Hernquist 
profile for the visible matter distribution, not the two component model used in the current 
fitting session. Nevertheless, one C3-I1 computG the mcitcliing distO/nccs for the obtained 
values of ri and rs for each galaxy, taking again a ~ ai ; the result is displayed in the last 
columns of Table |l]. 

Recalling the interpretation of Too as an indicator of the radius of the mimicked "dark" 
matter haloes. Table || clearly shows that the obtained n = —1/3 NEW haloes are always 
smaller than those conforming to the n = —1 isothermal sphere profile (except for NGC 
2434, which presents a best fit scenario with a single NEW "dark matter" component). 
Also, notice that r^o i is always larger than the range of the rotation curves, indicating 
that the isothermal sphere "dark matter" component extends further than observed; on 
the contrary, the value of Too 3 goes from about one half to less than twice the endpoint 
distance L of the rotation curves (roo 3/L = 0.57 for NGC 7145 and 1.73 for NGC 7507), 
indicating that the NEW halo is indeed less relevant than the isothermal sphere one in the 
composite scenario, where both contribute to the overall "dark matter" component. Even 
the odd case of NGC 2434, which displays no dynamically generated isothermal sphere 
component, yields a matching distance r^o 3 = 33.1 kpc of only about twice the endpoint 
of the corresponding rotation curve, L ~ 17 kpc. 

As expected, the matching with the cosmological background density for the n = —1/3 
NEW scenario occurs before that of the isothermal sphere, r^o 3 < 1, since the latter 
falls as r~^, while the former drops more steeply as r~^. Adding to this, the difference 
in scale between the characteristic densities pi » (since pi = 2KRi = 2k/ rf) should 
enhance this difference, not attenuate it. To ascertain this, one resorts to the dependence 
of the mimicked "dark" matter densities pdm 1 and pdm 3 on the latter: recalling Eq. 



(|4.12D , one has that 



1/2 1/2 fa\'^ 1/4 3/4 f'^\^ fr. (i\ 

Pdm 1 = Pi Pc' , Pdm 3= P3 Pc! (^-j , (5.6) 
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where pc = M/{2'ko?) = sets the visible matter density profile. Since pj^^ oc 1/ri and 
p]J^ cx one concludes that 
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- ^ 1/4 /a\ 1 a 



(5.7) 



Pdm 1 V«'^3/ ^rJ 10 r 

using Ts = 2GM and inserting typical values ~ 10~^ pc, a ~ 10 kpc, ri ~ 10 Gpc 
and r3 ~ 10^ Gpc. Hence, the NFW component of the mimicked "dark matter" typically 
becomes dominated by the isothermal sphere halo for r > a/10 ~ 1 kpc: this can be 
interpreted as an indication that the considered galactic rotation curves are "mostly flat" 
(again, with the exception of NGC 2434). 

5.3 Possible cosmological relevance 

The above discussion serves to assess the validity of the proposed mechanism of dynamical 
"dark matter" arising from a non-minimal coupling: the matching distances naturally 
assume values within the astrophysical context. This does not go against the previous 
claim that the obtained values for Ri and i?3 could give credence to the hypothesis of a 
relevant cosmological effect of the proposed model, thus enabling a possible unification of 
dark matter and dark energy. 

For simplicity, one retains the lengthscales ri and r-^ as fundamental, and convert the 
cosmological background density poo = >iH^ /SttG by defining the corresponding lengthscale 
(/9oo/2k)~^/^ = rn/V^, with rn = c/H = 4.2 Gpc the Hubble radius. With this notation. 



follows from the discussion of paragraph 4.3.1 that the non- minimal coupling Eq. (| 
reads, in a cosmological context, as 

A = + ^ f 'JL]"" U'JL\\ 10- + 10-', (5. 



J \n J 

inserting the typical values ri ~ 10 Gpc and ~ 10^ — 10^ Gpc. 

Given the form of the modification Eq. ( ^.l| ) , the above values hint that the n = — 1/3 
non-minimal coupling should not have a strong manifestation at a cosmological level; the 
strength of the n = — 1 power-law coupling, however, approaches unity and could lead to 
potentially interesting cosmological effects, i.e. a dynamically generated model for dark 
energy. Following the dominance of the gradient term in the r.h.s. of Eq. (|5.4| ) in the 
astrophysical setting, it should be remarked that both couplings could possibly be enhanced 
by the corresponding time-derivative terms in the Friedmann equations. This possibility 
will be developed elsewhere. 

5.4 Universality of the model 

As depicted above, the obtained values for ri and show a variation between galaxies 
that is clearly not compatible with the predicted universality: ri averages fi = 21.5 Gpc 
with a standard deviation ai = 10.0 Gpc, while presents = 1.69 x 10^ Gpc and 
(73 = 1.72 X 10^ Gpc. Although undesired, this is not unexpected, as there are several 
causes that might contribute to this deviation: although the selected type EO galaxies are 
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spherical at large, the small asymmetry due to the semi- axis difference 1 — 6/a<0.1 could 
translate into the difference in obtained values for rj. Furthermore, localized features and 
inhomogeneities of galaxies could contribute for the breaking of spherical symmetry and, if 
dense enough, serve as "seeds" for the dynamical generation of the mimicked "dark matter" 
component — or, due to the dominance of the gradient term in the r.h.s. of Eq. ( |5.4| ), if 
they vary significantly on a short lengthscale. 

Even if spherical symmetry were fully enforced, the real visible matter density profile 
might differ slightly from that obtained from the fitting of the visible component of the 
rotation curves, either due to a misconstruction of the latter or from an inadequacy of the 
choice of the Hernquist profile as a fitting function. Notice that this mismatch does not 
need to be large, so that the quality of the rotation curves reported in Ref. or of the 
related two-component Hernquist fits here used should not be held as an immediate culprit. 
However, one can envisage that a small deviation between the real and the fitted density 
profiles could include a large deviation of its first and second derivatives in very localized 
regions: again, since the r.h.s. of Eq. ( |5.4| ) is dominated by the gradient term, this could 
lead to a large deviation of the obtained "dark" matter component. 

Even if the simple modelling procedure followed in this study is sufficiently accurate, 
more fundamental issue could be responsible for the variation of the obtained values for ri 
and r^: first of all, the non-minimal coupling Eq. ( |5.lD could be too simplistic, and a more 
generalized model of the form 

f2{R) = \/5 [1 + 9i{R)] + § [1 + 92{R)] , (5.9) 

could yield the desired universality for the related parameters Ri and R^. In this sense, the 
current proposal aims at the the general depiction of the dark matter mimicking mechanism, 
and the functions gi{R) could be read from a more thorough adjustment with galactic 
rotation curves. 

One could also assume that the assumption of a linear curvature term fi{R) = 2kR 
might lead to the observed deviation from universality: indeed, power-law terms fi{R) oc 

as reported in Ref. [P, would also contribute to the dynamical generation of "dark 
matter". If the exponents are positive n > 0, this effect would be strongly felt at regions 
with higher curvature (i.e. density), closer to the galactic core: the failure to include these 
could perhaps account for the missing "dark matter" in the inner region of the obtained 
rotation curves. 

A second, perhaps more interesting issue could lie behind the obtained discrepancy: the 
non- minimal function f2{R) couples with the visible matter Lagrangian density, which is 
assumed to describe a perfect fluid form. Cm = —p- This can be relaxed in two ways: firstly, 
the perfect fluid model could be proven insufficient, and one should opt for a more adequate 
form Cm = Cmip) (or one depending on more variables related to the thermodynamical 
description of the matter distribution). 

Moreover, and perhaps more interestingly, one could invoke another related possibility, 
one with rather strong implications at a fundamental level: the non-minimal coupling does 
explicitly break the equivalence principle at large scales. This is not in conflict with labo- 
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ratory determination of this principle, which are eminently local and follow the reasoning 
that a violation of this principle is expected, for instance, due to the interaction of dark 
matter and dark energy [31, 32 1, or due to infrared gravity effects with some bearing on 
the cosmological constant problem [^]. Since this interaction (which might be modelled 
with non-trivial terms fi{R) or f2{R), as discussed before) should depend on the size of 
the galaxy itself, this could lead to the variation of the ri and rs parameters. 



6. Conclusions 



In this work one has addressed the issue of obtaining a solution to the dark matter puzzle, 
as embodied by the flattening of galaxy rotation curves. This is done by following the main 
phenomenological implications of models endowed with a non-minimal coupling of matter 
to curvature. 

In order to do so, one has first addressed the possibility of obtaining an adequate extra 
force from the non-conservation of the energy-momentum tensor of matter, for a completely 
flat and a decaying rotation velocity curve. In the former case, it was concluded that this 
required a logarithmic coupling of the form Xf2{R) = —v'^/mlog{R/R^:) (where m is the 
outer slope of the visible matter density p), which might by fitted by a simpler power- 
law Xf2{R) ~ (i?*/i?)°, with the asymptotic velocity given by = am. However, this 
solution suffers from a lack of flexibility that is not observed in nature, given the observed 
variety of the velocity dispersion of rotation curves, while the obtained solution implies an 
almost universal Voo- 

The second attempt was to solve this caveat, bearing in mind the phenomenological 
Tully-Fisher law, which roughly states a power-law relation between Voo and the total visible 
mass M: this is obtained by assuming that geodesical motion is indeed preserved, V^T^y = 
0, but that the metric itself is perturbed, thus providing the extra "force" required. In 
order to achieve that, one must solve the trace of the Einstein equations of motion, with the 
difference R — 2Kp providing the dynamically generated "dark matter" . We resorted to a 
power-law coupling with matter f2{R) = (R/Rq)^- since we expect the added dynamics to 
be manifest at large distances (vis a vis low density and curvature), this requires a negative 
index n. 

This approach yields a dark matter equivalent density directly related to a specific 
power of the visible matter density, thus accounting for the Tully-Fisher law in a natural 
way. Furthermore, the obtained dark matter component exhibits a negative pressure, a 
feature typical of dark energy models: this might hint at the potential of the considered 
model to unify both dark components of the Universe. 

Two different scenarios were considered: the mimicking of the NFW and the isothermal 
sphere dark matter profiles (n = —1/3 and n = — 1, respectively). Since a separate fit of 
the rotation curves with each power-law function f2{R) = (Ro/R)" did not produce results 
deemed satisfactory, a composite model for this non-minimal coupling was considered, 
allowing for a greatly improved adjustment. 

The value of the characteristic lengthscales ri and was obtained for each galaxy, 
instead of a simultaneous fit to all rotation curves: this yielded the order of magnitudes 
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ri ~ 10 Gpc and r^, ~ 10^ Gpc. This result allowed for the computation of the cos- 
mological background matching distances for each galaxy (which depend on their char- 
acteristic lengthscale) ; the obtained astrophysical range served to validate the proposed 
mechanism and conclude that the n = —1 isothermal sphere "dark matter" halo dominates 
the n = —1/3 NFW component. 

Furthermore, the n = — 1 scenario was also shown to possess possible cosmological 
implications, since ri ~ r//, the latter being the Hubble radius. This lends support to the 
proposed physical description, as well as to the hypothesis of a successful unification of the 
dark matter and dark energy components of the Universe. 

Given the lack of the desirable universality in the model parameters -Ri and i?3, several 
possible causes for the obtained variation were discussed: asides from the validity of the 
assumption of spherical symmetry, more fundamental issues were discussed — including 
a more complex form for the non-minimal coupling between geometry and matter or the 
curvature term fi{R), or for the Lagrangian density of the latter. A possible violation of 
the equivalence principle that depends on the physical size of a galaxy could also play a 
role in this variation. 

Although this work can be regarded as a first step, one concludes that the rich phe- 
nomenology stemming from a non-minimal coupling between matter and curvature enables 
an alternative description to dark matter. This should be contrasted with other candidate 
models, namely the MOdified Newtonian Dynamics (MOND), which is by itself purely 
phenomenological, or based upon an extensive paraphernalia of vector and scalar fields 



(the underlying Tensor- Vector-Scalar (TeVeS) model) [^, 35, 36, ^ (see Ref. [38] for a 
critical assessment). 

Finally, we remark that the proposed model should not be held as a competing pro- 
posal to standard f{R) theories (that is, with fi{R) = f{R) and f2{R) = 1), but as 
complementary: it has been shown that a general action with non-trivial fi{R) and f2{R) 
dependences is equivalent to a scalar-tensor theory exhibiting the proposed non-minimal 
coupling (as studied in Ref. [^] and discussed in Appendix B), so that both may be re- 
garded as limiting cases. Indeed, in the metric formalism (where the Christoffel symbols 
are assumed to depend on the metric), f{R) gravity is equivalent to a Brans-Dicke (BD) 
theory with a vanishing BD parameter a; = and a suitable potential, together with a 
non-minimal gravitational coupling of the BD scalar field. 

From a physical standpoint, the possibility of describing galaxy rotation curves via a 
non- minimal gravitational coupling (as obtained here and in Refs.|8|, P) may be interpreted 
straightforwardly: traditionally, the discrepancy between the observed velocity dispersion 
and the lower predicted value is due to the addition of some missing "dark" matter form. 
By resorting to a non-minimal gravitational coupling, either explicitly or via the discussed 
scalar-tensor equivalence, one is instead asserting that this disagreement is solved not 
by an additive mechanism, but a multiplicative one: the influence of the existing visible 
matter is suitably increased in the outer galactic region by the non-minimal coupling. The 
identification of this heightened effect as dark matter, albeit artificial, then serves the 
purpose of contextualizing the obtained results within the vast literature on the subject. 
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A. Power-law density profiles 



A plethora of different density profiles exist for both visible as well as dark matter com- 
ponents, arising both from observations as well as N-body simulations. The key feature of 
the proposed approach is the ability to mimic dark matter in the outer region, where the 
curvature is low enough so that the effect of the non-trivial, inverse power-law gravitational 
coupling f2{R) becomes manifest. For this reason, it is not necessary to fully account the 
specificity of each of the competing models: it is sufficient that these exhibit a behaviour 
which is dominated by a power-law for distances above a certain threshold a. 

Several models exhibit this feature: amongst these, one of the most discussed are the 



so-called generalized (spherical) cusped profile [p^, HO, 41, 42], given by 



where pcp sets the density scale, a the length and 7, m are the inner and outer slopes, 
respectively. 

Typical cusped density profiles include the 7 = 1, m = 3 NFW profile for dark matter 
j43U , and the 7 = 1, m = 4 Hernquist profile for visible matter Q]. The latter yields a 
finite total mass, M = 27ra^pcp; substituting into Eq. ( |A.lJ ), one also gets the mass enclosed 
within a sphere or radius r, 
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so that half of the total mass M is contained within /2 = (1 + V^)a. 



1/2 

A more general model is the so-called (a,/?, 7) model |^^, where the density is given 



by 



, r 

l+(- 

a 



(7-/3)/a 



(A.3) 



which, for r ^ a, assumes the simpler power-law p(r) oc . 

The above models are all divergent as r — )• 0, thus to alleviate this and accommodate 
the apparent fiat density profile observed near the core, some other proposals have been 



put forward, including the Burkert profile |22], given by 



Pir) = Pv- 



(A.4) 



' (r + a)(r2 + a^) ' 

which, in the outer region r ^ a, behaves as p{r) oc r~^. 

All the above models present a shift in behaviour at r ~ a, changing the slope of the 
inverse power-law from (assuming the notation of the Hernquist profile) 7 to m. It should 
be remarked that the mechanism considered in this work for mimicking the outer density 
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profile of dark matter is only valid for this region, since only by accident can the scaling law 
Eq. ( 4.121) yield the r"™"' behaviour for the dark matter component from the r"™ power 



law for visible matter, and simultaneously produce the law for the former, from the 
r^'^ profile for the latter (this is the case if one assumes the Hernquist profile for visible 
matter and the NFW profile for dark matter, as both exhibit 7 = 7' = !). This caveat from 
the current proposal, although relevant, does not hinder its main advantage: to reproduce 
the dark matter halo in the outer region, where it is most relevant. 

Finally, one presents what is arguably the simplest of dark matter density profiles, the 
singular isothermal sphere, given by 

where foo is the asymptotic velocity in the flattened region ( i.e. the velocity dispersion). 
Indeed, this allows for a straightforward computation of a flattened rotation curve, with 
v(r) = Voo everywhere. 

The isothermal sphere may be regarded as a generalized cusped profile Eq. ( [A.lj ), 
with an inner slope 7 = 0, outer slope m = 2, and vanishingly small a — t- 0, while keeping 
Pcpa^ = vI/{27tG). 

One might consider a different family of density profile models, which do not assume 
an inverse power-law behaviour at all, but instead rely on a modified exponential decay, 
usually depicted through one of two closely related models: the Einasto R^^p density profile 

lie 



p[r) = exp <^ -dp 



a/ 



(A.6) 



where dp is an adequate function of p, such that r = a encloses half of the total mass; and 



the deprojected Sersic law |^ 



Pi^') = Pv ^ exp 



dm I 



7- \ 1/m 

a 



(A.7) 



with and d^ functions of m. 

Most simulations assume a power-law density for visible matter (most often, the Hern- 
quist profile), while the Einasto and deprojected Sersic law compete with the NFW density 
profile in describing dark matter halos. Simulations using a deprojected Sersic profile pre- 
dict that the exponents for visible and dark matter are = 4.29 and mdm = 2.96, 



respectively [|50 |. 



It is straightforward to see that exponential density profiles such as those depicted 
above would also be obtained for the dark matter component, if one assumed similar model 
for visible matter to begin with; however, the scaling law Eq. ( [4.12| ) will only change the 
overall value oi pm and dm, through 

Pdm=P*r-^ — , ddm = d^:—^ — , (A.8) 

1 — n 1 — n 
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while the exponent 1/m will remain unaltered by the scaling. Given that pm and dm are 
directly dependent on this quantity, one concludes that the proposed mechanism cannot 
satisfactorily reproduce the above exponential density profile for dark matter starting from 
a similar model for visible matter. This is a caveat of the current work — although it 
should be stated that the specific choice f2{R) = {R/Rq)^ was tailored to deal with the 
power-law scaling, and a more appropriate choice might resolve this issue. This will not 
be discussed here, since the underlying idea would remain the same: to mimic dark matter 
through an appropriately driven dynamical effect arising from the non-trivial gravitational 
coupling. 



B. Equivalence with multi-scalar theory 

Following Ref. p|, one may rewrite the two couplings obtained in sections ^ and Q as a 



multi-scalar theory |51] with action 



S 



I 



-gd X 



where (p^ and 93^ are scalar fields, aij is the field- metric 



(B.l) 



0"< 



1 1 

-1 
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the potential is given by 



exp 



2^/3 1 



2k 



exp 



2V3 1 
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and £.* = exp 

The two scalar fields are related with the scalar curvature and the non-trivial fi{R) 
and f2{R) functions: 



1 

f = ^ log 



FiiR)+F2iR)C 
2k 



ip' = R, 



(B.4) 



where Fi = dfi/dR. Inserting C = p, fi{R) = 2kR and f2{R) = {R/RqY, one gets 



V5 



V3, 
— log 

2 ^ 



\RoJ Po 



^log [l + ng-^] 



(B.5) 



Hence, one concludes that the dimensionless function g not only simplifies treatment of Eq. 
( |4.6| ) , but embodies the added degree of freedom resulting from considering a non- minimal 
coupling between matter and curvature. 
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Figure 1: Left: Observed rotation curve (daslied full), decomposed into visible (dotted) and 
dark matter (dashed grey) contributions |27| , superimposed with the mimicked dark matter profile 
(full grey) arising from the composite non-minimal coupling and resulting full rotation curve (full). 
Right: Log-Log profile of the visible matter density (dotted), mimicked "dark matter" contribution 
(dashed) and the sum of both components (full) . 
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Figure 2: Continuation of Fig. 1. 
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Figure 3: Profiles for visible matter density (dashed), total density (full) and contributions arising 
from the gradient (full grey) and static (dotted) solutions with n = — 1/3 (left) and n = — 1 (right), 
for NGC 7145. 
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Figure 4: Observed rotation curve (dashed full) , decomposed into visible (dotted) and dark matter 
(dashed grey) contributions |2^, superimposed with the mimicked dark matter profile (full grey) 
arising from the non- minimal couphng /2(i?) = {/R3/R (left) and /2(i?) = Ri/R (right). 



-39- 




10 20 30 40 50 10 20 30 40 50 



Figure 5: Continuation of Fig. 4. 
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